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Abstract. We present results on both the intensity and phase-dynamics of the 
transient non- linear optical response of a single quantum dot (SQD). The time 
evolution of the Four Wave Mixing (FWM) signal on a subpicosecond time scale is 
dominated by biexciton effects. In particular, for the cross-polarized excitation case 
a biexciton bound state is found. In this latter case, mean-field results are shown to 
give a poor description of the non-linear optical signal at small times. By properly 
treating exciton-exciton effects in a SQD, coherent oscillations in the FWM signal are 
clearly demonstrated. These oscillations, with a period corresponding to the inverse of 
the biexciton binding energy, are correlated with the phase dynamics of the system's 
polarization giving clear signatures of non-Markovian effects in the ultrafast regime. 
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1. Introduction 

Exciton dynamics experiments are attracting continuos interest because of their 
suitability to explore dephasing effects of single exciton and multiexciton complexes [UEj- 
Furthermore, recent proposals for solid state quantum computing systems [3J EJ G3 Ej, 
have stressed the importance of properly controlling multiexciton coherences. Recently, 
the non-classical behaviour of the light emitted from a SQD, photon antibunching in 
the fluorescence spectrum, has been observed in the ultrafast regime jH]. Additionally, 
the search for single photon sources jH| has triggered experimental interest in transient 
multiexciton coherences. All of these exciton based phenomena reflect the importance of 
interparticle interactions (electrons and holes) and their couplings to the environment, 
e.g, phonons. Both coherent and incoherent effects are experimentally detected. 
Incoherent effects are related to the system-bath interaction and from a theoretical point 
of view they are usually modeled within a Markov approximation. Coherent effects are 
associated with multiple particle correlations and can be observed in the subpicosecond 
time scale. Consequently, a detailed understanding of the coherence persistence at short 
times is of most importance in nanostructure systems. 

Time resolved FWM experiments give signatures of particle-particle correlations 
and particle-environment interactions. In a typical FWM experiment, two simultaneous 
excitation pulses co- or cross-circularly polarized, with wave vectors k\ and k 2 , propagate 
onto the sample. A third laser pulse, with wave vector k^, is sent time delayed at 
t 2 = T. For moderate excitation intensities, exciton-exciton interaction arises. The 
time evolution of the third order polarization P^ 3 \t,T) = \P^ 3 \t, T^e 1 *^' 7 ^ describes 
the properties of the diffracted light in the k^+^ — ki direction. At very short times after 
the excitation, particle correlations build up. These correlations have been observed 
using spectral interferometry techniques in bulk systems [TU]. It is also possible to 
study the polarization phase dynamics itself, $(t,T). Moreover, other experiments for 
measuring p( 3 \ like pump and probe schemes, have also been performed. In particular 
in Ref . 0] a simple four level model has been used to get the non-linear optical response 
(x^) °f a SQD. From the frequency spectrum of x i t ne possibility of getting exciton 
entanglement in a SQD has been predicted. However, in that work exciton correlations 
were included within a Mean Field Approach (MFA). 

Exciton dephasing times (T 2 ) have been well characterized in bulk and quantum 
well systems. In particular, ultrafast dynamics experiments, like time-resolved speckle 
analysis, have been used to get dephasing times information in solids [HI [12J [T3j. On 
the other hand, recent experiments in bulk systems show the build up of many particle 
correlations in the femtosecond scale of time|14j. Recently, all of these dynamical 
phenomena have been observed for first time in Self Assembled QDs using transient 
FWM spectroscopy where the time integrated FWM signal oscillates with the 
biexciton binding energy. Besides that, dephasing mechanism from Optical Longitudinal 
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(LO) phonons has been observed in CdSe Bulk and quantum dots systems, showing 
clearly non-markovian effects associated to the coupling between carriers and LO 
phonons [Iff], Due to its atomic like characteristics, a SQD offers the unique possibility 
of manipulating the number of particles and their Coulomb interactions (exciton-exciton 
(X-X)) with direct consequences on decoherence processes control. Thus, rather long 
dephasing times (T 2 ~ 40 ps) have been reported fTf\ H%1 IT§]. 

A natural question arises when ultrafast scale of times are involved: what kind of 
quasiparticle (free electron-hole pairs, exciton or bound excitons) are dominant?. In 
order to answer this question we consider the Phase Space Filling (PSF) effect and X-X 
interactions on the same footing. Within MFA, the FWM signal shows contributions 
coming from both, PSF effects as well as renormalized X-X interactions. However, 
this approximation fails to explain correlations at very short times. Therefore, MFA 
should be improved by considering quantum fluctuation effects to explain transient 
FWM experiments. In particular, more elaborated theoretical approaches like the one 
reported in Ref. |2()j. where a truncation in the hierarchy of polarization equations of 
motion to fifth order in the optical electric field is proposed, are shown to improve 
MFA results by including four particle correlations. Aditionally, some previous reported 
works [2QJ |^] show that in bulk systems a mean- field approach to calculate the FWM 
signal gives errors because the exact X-X correlations are neglected. 

In the present work, we go beyond previous limitations like few level systems and 
MFA. We use a truncation scheme similar to one developed for bulk systems[2J, in 
which the contribution from exciton and exciton-exciton effects can be split and treated 
on an equal footing to any desired order in the optical field [2D]- However, the truncation 
scheme fails to evaluate exactly the X-X Coulomb interaction due to the intensive 
computational work involved in its application to bulk systems. In order to get feasible 
results, this scheme has been used mapping the problem to a one-dimensional Hubbard 
model [21]. In the present work we undertake this problem and we include exact X-X 
correlations in order to find the non-linear optical signal, intensity and phase dynamics 
in a realistic model of a SQD in the ultrafast regime. In Section 2, we briefly review 
the theoretical background on which our results are based. Section 3 describes the main 
results of this work. Our conclusions are summarized in Section 4. 

2. Theoretical Model 

Our starting point is the Hamiltonian for a system of electrons and holes in a two- 
dimensional parabolic quantum dot excited by external laser pulses 

H = V E e c f c + V E h h ] h 

+ yi »t£(tHA,s< + E it:Myh v ,s>c v ,s 

v,rj(s,s') u,a(s,s') 

~ £ Y(^i,^2\V e - h \u 3 ,^)ci 1>s hl 2>s ,h V37S 'C V4tS 

s s ' 

"3 1^4 
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Z "1>"2> s,s' 

+ o E £("l> ^l^fc-hl^ ^) h l^,s h ^,s' h ^,s'h VA ,s (1) 

"3 '"4 

where c(/i)J, s , c{h) VtS creates and destroys one electron (hole) in the state v (labeled by 
quantum numbers (n e (/0> m e(h))) with spin s (f , spin up, j, spin-down). V^_ e , V^_/j, 
and V^-h, denote electron-electron, electron- hole and hole- hole Coulomb interactions, 
respectively. Single particle energies for electrons (holes) are given by = (n e (fe) + 
m e (h) + l)^e(fe), with confinement energies u e (h) = l/ m «)C(k) = 1)> where l e ( h ) is 
the parabolic confinement length size of the electron (hole). The set of laser pulses used 
to excite the SQD is described by the envelope amplitude £{t) and associated dipole 
moments which take into account pulse polarizations. The spins of the electron 
(s) and hole (s') determine the corresponding polarization index a i.e. a = a(s,s'). 
By diagonalizing the time-independent part of the Hamiltonian Eq.Q (S(t) = 0), the 
energies and wave functions for one and two excitons in a SQD are obtained, from which 
the non-linear optical response is calculated. 

Exciton and Biexciton wave functions are given respectively, in terms of non- 
interacting electron-hole pairs, as 

= E K«(i,X,J°> ( 2 ) 

"3 1 "4 

|0) is the SQD ground state and U2 and ^ 4 are the exciton and 

biexciton amplitudes, respectively, which take properly into account the carriers spin 
configuration. In the following formulas we will stick to the notation of Eq.(J2J) where 
the index n labels an exciton state, whereas the index m is reserved for the labeling of 
biexciton states. 

In discussing the time resolved spectra, we will be concerned with the situation 
in which the SQD is excited by very short laser pulses (5 functions) in resonance with 
the Is heavy-exciton state. Therefore, a big number of single particle states has to 
be included in the calculation of one and two exciton spectra. The applicability of 
the 5-pulses limit can be extended to typical accoustic periods when interaction with 
LO phonons can be neglected, compared with the deformation potential. However, it 
will be very important in II- VI semiconductors due to this high enough polar coupling. 
Our results are related with III-V systems, in which this coupling is low. In this way, 
when coherent phenomena coming from the oscillations of exciton or biexciton-phonon 
complexes are lost, incoherent phenomena aquires more importance and after this time 
it is well known that a phenomenological model of decoherence, as the one we will adopt 
here, can be used. 

We perform our theoretical calculations following a truncation scheme of the 
Hilbert space of one exciton and two interacting excitons. After a long algebra it is 



Ultrafast non-linear optical signal from a single quantum dot: exciton and biexciton effects^ 



possible to give analytical expressions for the transient non-linear optical signal, which 
is proportional to P^ . It can be described as the sum of three terms: one representing 
the PSF (P PSF ) contribution, another one corresponding to the MFA term (P MF ) 
and a third one related to the X-X contribution {P F )- This last term describes exact 
four-particle correlations and memory effects in terms of the so called force-force time 
correlation function. All these terms can be calculated for a SQD described by the 
Hamiltonian Eq. JIJ , giving place to expressions similar but not equal to those of Ref . [21] 
(the differences are due to the confinement geometry of the dot, see the discussion in 
the next section). The resulting expressions for each one of these contributions for a 
SQD are given explicitly by 



rT 0,< T 2, rT 3 



where a n>(7 is related to the exciton wave function at zero relative distance, '^o-nl'al 
is the PSF parameter (for notation details, see Ref.j2I]) and 0(x) denotes the step 
function (O(0) = l). The MFA contribution can be written as 



n ,n 2 ,n 3 , ^x,no ' ^x,n± ^x,n 2 ^x,n 3 > ^v*- 



o "0. o "2. CT 3 



where 



and 



x [6(-T)e(t + T)A x (t, T) + e{T)Q(t)A x {t, 0)] (4) 
= E0"«,m - <**,r« - "*,ni)B2%\{rn)B%%l{m) (5) 



A (t T^j = < ^~' r ^ LUx ' n ^~ UJ x,n 1 ^^1,112 ^^1,113 ))* g(^ — * \^x,tiq +Wa!,ni — ^x,n 2 ~ ^x,n 3 JjT ^Q 'j 

Similarly, the X-X contribution is given by 

(t, T)e(— +r )' = ia ni , ai ( £ P ni (m)) - P^ (t, T)e^^ (7) 



in 



where 



p ( m ) = a a* a* B r7 °' ai (m)B CT2 > a3 (m)- 1 ,m — 

^niV rl ) a n ,cT ^n 2 ,a 2 ^n 3 ,a 3 n n Q ,n 1 \' n ) ID n 2 ,n 3 \l n ) _ _ _ T 

HQ, 112,13, LUxx,m LUx,no '>'- xx 



r7 0< a 2-< 7 3 

x [(G(-T)e(t + T)A xx (t, T) + Q(T)Q(t)A xx (t, 



with 



A (t X 1 ) — ^ xx ~' r ^ LUx ' n ^^x.ni — ^xx,m))t xx ~ i(^x ,n +^x ~ ^>xx ,m ))T 

In Eqs.(J5J) and (J7J), m runs over all (bound and unbound) biexciton states. Bound 
biexcitonic states are those for which its energy lies below two times the ground state 
energy of an exciton. Aditionally, P^'Z^nl'vl an d Pn 1 ,ai( , m) take into account the 
exciton and biexciton transition weights on the signal produced by the n\ exciton 
state. The exciton transition weights Bnlln^m) are defined to be the matrix elements 
(0\B nu(Ti B n v \m) xx of the product of two n th exciton destruction operators B na 
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between the exciton ground state (0| and a biexciton state \m) xx . Finally, all we need 
is: exciton, uj x>n , and biexciton, uj xx>m , energies and their respective wave functions, 
which can be numerically obtained to any desired precision in the SQD case. From 
these ingredients, the exciton transition weights can be calculated as _B„ o ni (m) = 
(0\B no B ni \m) xx where the exciton destruction operator is B n ^ a = ^Z Vl V3 U2 c Vl)S h V2 ^i 
(a = a(s,s')) and is the n-th exciton amplitude (cf. Eq.(J2J)). One of the main 

results of our paper is that, in contrast to an extended system, for QDs the term 
Pnl^ Tl does not vanish, even in the case of cross-polarized (CP) excitation. If the 
dot confinement energies are reduced, in the CP case ^ decreases in magnitude, 
approaching zero in the limit of zero confinement (extended system). This can be seen 
by direct evaluation of the parameter j3 in Eq.fjSJ), that gives information on the strength 
of the signal as predicted by MFA. Although, the (3 term contains information about 
the renormalized biexciton energy, this term is not responsible for oscillations and only 
contributes with an spectral weight given by the Bnlln] changing the intensity of the 
FWM signal. 

The total dynamics of excitons and biexcitons comprises the MF and the exact X-X 
contributions. The temporal dynamics dependence of P MF (apart from an exponential 
damping) is contained in Eq.© where only contributions from differences of exciton 
energy levels appear, which are negligible in the case of resonant excitation. On the 
other hand, the temporal dependence of the X-X contribution is dominated by energy 
differences between two exciton states and one biexciton state, e.g. Eq.Q. In particular, 
in the CP bound biexcitonic state is present, causing oscillations in the signal 

with a frequency corresponding to the energy of the bound state. 

Exciton and biexciton dephasing rates associated to non-radiative mechanisms are 
described (phenomenologically) by T and T xx , respectively. In the following we will 
concentrate only in the lowest exciton state, i.e. ri\ = n 2 = n 3 = = Is with 
T xx = 2T. In our model, the dephasing time T 2 is defined as the inverse of T xx . 

3. Results and discussion 

In order to calculate the FWM signal for a typical self-assambled SQD 22, 23 , we 
consider a Jno.5Gao.5As quantum dot with e = 12.5, m* = 0.068, m* h = 0.2, oj e = 20 
meV and u>h = 5 meV. To guarantee a good convergence for the spectrum of zero total 
angular momentum excitons and biexcitons, we diagonalize a 72x72 Hamiltonian matrix 
for excitons and a 3250x3250 Hamiltonian matrix for biexcitons. The exciton spin is 
determined by the circular polarization of the excitation pulses. Further in this work 
we discuss only the cross polarization excitation case. With this set of parameters the 
biexciton binding energy is -1.27 meV 

As discussed in Section 2, the FWM signal clearly depends on T, on the exciton and 
biexciton energies as well as on the transition spectral weights. Excitation pulses create 
coherent excitons whose dephasing time (T 2 ) can be monitored by the delayed pulse. 
As the time delay increases, the FWM signal decreases from which T 2 can be obtained. 
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With these ingredients, and considering that more than one electronic excitation can be 
induced by the laser and four particle correlations occur for times lesser than T 2 /4, we 
expect that interference phenomena appear in the ultrafast time scale. Therefore, the 
coherent FWM signal and the phase dynamics give clear signatures of exciton-exciton 
interactions in the ultrafast regime. In the following we will not be concerned with T2 
determination, which is to be taken as a free parameter. Pulses with zero time delay, 
T = 0, will be considered. 

In the weak intensity regime both PSF and X-X effects determine the nonlinear 
optical properties. We will concentrate on these two contributions. PSF contributions 
to the FWM signal are only important in the parallel spin exciton case, where X-X 
effects are negligible, because a bound biexciton state does not exist at all. By contrast, 
for the cross polarization situation we are interested in, PSF effects are vanishing while 
X-X effects become dominant. It is important to remark that in extended and isotropic 
systems the MFA gives a zero contribution to the FWM signal, because only excited 
excitons with zero center of mass momentum contribute, i.e.,(<f = 0) |21j. However, in 
SQD this condition can be relaxed due to the fact that the parabolic potential absorbs 
the incident momentum and excitons with center of mass momentum different from 
zero can be excited. This point becomes particularly clear by considering the mean 
field parameter (3 defined in Eq.©. Whereas for extended systems this parameter is 
exactly zero [21], for a SQD it becomes a function of the confinement energies 0J e {h). 
The dependence of (3 on u e is depicted in Fig. (1). We consider a constant relationship 
between electron and hole confinement lengths, i.e, l e = lh- It can be seen that in 
the limit of small confinement, i.e. pure two-dimensional systems (o; e> ^=0), (3 tends to 
zero, as it should be. This means that the usual condition of (q 7^ 0) for the carriers 
Coulomb interaction in extended systems is broken in confined systems. We have tested 
our results taking a huge set of excitonic states and we found that higher energy states 
give a negligible contribution to (3. 

The phase-dynamics and the non-linear optical signal are calculated for different 
decay rates. In order to better clarify, why MFA is not adequate to describe non-linear 
optical signals at small times, we plot in Fig. 2, separately, the amount of particle 
correlation contribution to the FWM signal as described within MFA (EqJH) and exact 
X-X effects (EqEJ), which include memory (non-Markovian) effects, for T xx = 0.125 
meV. The usual and simplest approximation for this last term yields to the Markov 
approximation and MFA results. By contrast, in this work we have performed a full 
numerical evaluation of X-X effects by summing over all biexciton states. Hence, we 
plot the absolute values of both, the real part of P*^(t, T)e (i< ^"i +r) ' (MFA term) and 
the real part of P£ ai (t, T)e( laJ:c n i +r ^ (X-X term). The MFA result grows monotonically 
with t, whereas the X-X result shows beating oscillations with a period corresponding 
to T xx = 1ixoj~ x x fiVV « 5ps, the biexciton binding energy. As can be seen, non-Markovian 
effects arising from exciton-exciton correlations are the main source of discrepancies 
between MFA results and the exact ones. This means that ultrafast spectroscopy in 
SQDs should register memory effects associated to four particle correlations. 
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In the inset of Fig. 2, we show the spectral weight B(m) as a function of biexciton 
energy levels (the zero of energy corresponds to twice the single ls-heavy-exciton 
energy). We have only plotted those levels that have total angular momentum L z = 0. 
The bound biexciton level is clearly seen at negative energies. Although, unbound 
biexciton states have non-negligible spectral weights the bound biexciton still has a 
comparable weight. As can be seen in Fig. 2 MFA fails to describe properly four particle 
interactions at very short times, where particle-particle correlations are important. This 
term does not oscillates and only average the exact result. Besides that, the contribution 
in Eq|?l takes into account a higher degree of correlation going beyond MF and gives 
oscillations at very short times. This oscillations provide information about strong 
correlations between excitons. Thus, the bound biexciton dynamics is directly linked 
to non-Markovian effects. By contrast, in the long time limit (t > T2), memory effects 
are lost and MFA results approach to the real part of the first term in EqEI The FWM 
signal becomes dominated by non-radiative dephasing effects, producing in this way a 
simple exponential decay. 

In order to assess how the main features described above are reflected on the optical 
dynamical characteristics, we plot the non- linear FWM intensity |P (3) (t, 0)| in Figs. 3(a) 
and 3(c), and the dynamical phase $(t, 0) in Figs. 3(b) and 3(d), for dephasing rates, 
T xx = T 2 ~ =0.125 meV and 0.5 meV, respectively. To better illustrate the results, 
a comparison is made between MFA and X-X results on the same plots. For a low 
dephasing rate (T 2 ~ 8 ps), strong oscillations are clearly seen in Figs. 3(a) and (b) 
for t < T 2 . The X-X interaction, induces a 7r phase shift with respect to the incident 
electric field. The phase of the signal oscillates during a typical time T 2 . FWM 
signal shows a strong peak at t « T 2 /4 whereas it should have a maximum at t ~ T 2 
from MFA. For high dephasing rates (T 2 w 2 ps), Figs. 3(c) and (d), the attenuation of 
both the intensity and phase of the third order optical signal is evident. At this time 
scale, the T 2 dephasing time is shorter than the time for which coherent effects could 
be seen. It is worthnoting that while the phase shows a similar starting behaviour as 
compared with the low dephasing case, the FWM signal intensity is drastically reduced 
by roughly a factor of five (notice the change of vertical scale). 

Most importantly, the phase dynamics exhibits novel features at small times when 
exciton-exciton correlations are considered. By contrast, MFA results yield to a real P^ 
for any time. Therefore, no phase dynamics is observed within MFA. In particular, X-X 
effects present a correlation between the first FWM signal minimum and a ir jump in its 
phase. For long times the phase goes to a constant value. However, after the extinction 
of the FWM intensity, there is no longer a clear physical meaning for this phase. Hence, 
oscillations in the phase dynamics, with a beating frequency controlled by the biexciton 
binding energy, could bring important information about exciton correlations in SQDs. 
In contrast with higher dimensional systems (quantum wells and bulk samples), a SQD 
provides us with an adequate system where exciton and biexciton discrete energy spectra 
can be tailored by changing the confinement potential. In this way, optimal conditions 
to enhance phase memory could be explored. 
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4. Conclusions 

In summary, we have studied the non-linear optical response from multiexciton 
complexes in a SQD. For cross-polarized excitation, the creation of bound biexcitons 
is possible. They dominate the non-linear polarization dynamics in the low density 
regime. The phase dynamics and the FWM signal in the ultrafast scale of time have 
been obtained. By including all the numerically determined biexciton states, exact 
exciton-exciton correlations are evaluated allowing to go beyond simple MFA results. 
In particular we found that the phase dynamics and the FWM intensity oscillate with a 
period which coincides with the inverse of bound biexciton frequency. We demonstrated 
clearly that the usual Mean Field theory, which is closely related to the Markovian 
approximation, does not describe correctly exciton correlations at very small times. 
However, our results show that at long times, MFA results for both intensity and phase 
dynamics, are close to those produced by a more realistic calculation. 
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Figure Captions 

Fig. 1: MFA contribution {(3 term) as a function of the electron and hole confinement 
frequencies, for T xx =0.125 meV. 

Fig. 2: Absolute values of the MFA contribution (real part of P^ F (t, T)e ( ^-"i +r) *) and 
X-X contribution (P£(t, T)e (i ^-"i +r) ') to the FWM signal, for T xx =0.125 meV. Inset: 
Spectral weight of the biexciton wave function. Bound biexciton energy =-1.27 meV. 
The long-dashed line denotes the zero energy. 

Fig. 3: 1^=0.125 meV: (a) Time resolved FWM signal; (b) Phase dynamics. T xx =0.5 
meV: (c) Time resolved FWM signal; (d) Phase dynamics. 



